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Abstract 

We analyse a model where the anomalies of the bond-stretching LO phonon mode are caused 
by the coupling to electron dynamic response in the form of a damped oscillator and explore the 
possibility to reconstruct the spectrum of the latter from the phonon measurements. Preliminary 
estimates point to its location in the mid infrared region and we show how the required additional 
information can be extracted from the oxygen isotope effect on the phonon spectrum. The model 
predicts a significant measurable deviation from the "standard value" of the isotope effect even if 
the phonon frequency is far below the electron excitation spectrum, provided the latter is strongly 
incoherent. In this regime, which corresponds to the "mid infrared scenario", the phonon linewidth 
becomes a sensitive and informative probe of the isotope effect. 

PACS numbers: 74.72.-h, 71.38.-k, 73.20.Mf 
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I. INTRODUCTION 



One of the reasons the lattice degrees of freedom of cuprate superconductors remain 
an intensively debated topic on both experimental and theoretical sides is the existence of 
strong anomalies in some of the phonon spectrum branches and their possible relation to 
electronic properties Q, [3, [s, |4|. In particular, the high frequency, about 80 meV, in-plane 
bond-stretching LO ( BS or the "half-breathing") mode shows large damping, typically sev- 
eral meV, and dispersion dip of about 20% for momenta (qy) around 0.25 r.l.u. The 
line broadening and softening increase with hole doping and are sensitive to temperature 
close to T„ (see jst for a review; the superconductivity induced anomalies have been dis- 
cussed in |6[ for Raman active phonons). The first measurements of the Cu-0 BS mode 
only reported a cosine-like dispersion, which can be understood in terms of conventional 
calculations. Subsequent measurements, however, pointed to anomalous behavior related 
to charge inhomogeneities and different interpretations have been given. McQueeney and 
co-workers Q reported a dynamic doubling of the unit cell in the Cu — O bond direction 
and discontinuous phonon dispersion. This has been related to charge stripes with a slower 
dynamics than the BS phonon frequency and formation of a second phonon mode with a 
lower frequency. In contrast, Pintschovius and Braden have found a smooth BS mode 
dispersion with an anomalously large slope and broadening around q = (0.25,0,0). Reznik 
and co-workers [9i] have studied La2-xBaxCuOi for x ~ 1/8, where quasistatic stripes ex- 
ist, and have assigned the observed strong enhancement of the phonon broadening to the 
sharp softening. D'Astuto et al. \lO\ have carried out high precision X-ray scattering ex- 
periments on La2-xSrxCu04+s and concluded that this phonon broadening is caused by 
intrinsic damping and not to an apparent one due to steep dispersion. 

In an earlier paper based on a semi-phenomenological model, we have suggested 
that the anomalously strong softening and broadening of the half-breathing mode could be 
related to the other anomaly observed in the mid-infrared region of the optical conductiv- 
ity spectra (MIR) in high-Tc cuprates 13]. Apart from making this connection, the way 
our model addresses the problem of the phonon anomaly is opposite to other approaches: 
phonon properties are taken as an experimental input which then serves to determine cer- 
tain properties of the electron excitations. It explains that the measurements taken so far 
are insufficient and motivates new experiments. Namely, to bridge the existing gap it is 
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important to know the effect of oxygen isotope substitution upon the phonon spectrum. 
In what follows we describe the details of this approach which can be directly applied to 
the analysis of such experiments as well as to test the consistency of the model. We also 
calculate the expected values of the measurable parameters depending on the frequency of 
the electronic oscillator which spans the range between MIR and near resonance frequen- 
cies. As the BS phonon mode couples to electron density excitation only at g 7^ 0, it could 
give a unique opportunity to access the finite momentum evolution of the MIR excitation. 
On the other side, the present work could stimulate further theoretical investigation of the 
isotope effect within microscopic models, as for instance on the mechanism of dynamical 



stripe fluctuations 13(]. An interesting possibility is implied by the analysis of site-selective 



isotope substitution on Raman active phonons 



14l | when, in contrast to total substitution. 



a significant change of coupling occurs due to modification of phonon eigenvectors. 
II. THE MODEL 

Description of the phonon softening relies on the Dyson equation for Green's function 
(GF): ^ 

where A is the dimensionless coupling constant and P {q, u) electron density response func- 
tion [L5\. It is assumed that the part of this response which is relevant for the phonon 



softening can be considered in the form of a damped oscillator 11 1 



The oscillator is parametrised by the full- width at half-maximum (FWHM) Tq, frequency 
Qq and oscillator strength rjq which have to be determined by analyzing the data on the 
anomalous BS phonon branch. The experimentally determined values of the linewidth 7^ 
and frequency ujq are used to define the standard form of the renormalized phonon GF 

In our model the phonon GF is found by solving for the complex poles of u = Uq — ijq : 

2 2,2 ^9^0 K - r,7g - + 7') 

^0-^g+lg = -~ 



{nl-Tq^q-U;^q+J^y+i{Tq-2^q)u;qy 
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(r, - 27, 



(3) 



Thus, we have the correspondence ujn 



47g and ^qUJq = 2co'g7g, where the phonon 



hnewidth measured in experiments, 7^, is roughly twice as large as the quantity we use in the 
calculations, 7^. Still these two equations are insufficient to determine the three parameters 
of the model, Tg,flg and = Ar/g sin^ (g^), and we further discuss additional input that 
could be obtained by varying the frequency of the phonon mode, i.e., by oxygen isotope 
substitution O^^ O^^. 

In fact, there exists a detailed information about the isotope effect (IE) on electronic 

n □ 

properties of cuprates, e.g. on Tc [16[ and ARPES [17j, but there have apparently been no 
attempts to consider the effect such substitution has on the phonon spectra. It is known that 
the high energy phonon modes, and in particular the one considered here, are dominated 
by oscillations of oxygen ions. Therefore such experiments should contain important infor- 
mation about both lattice and electronic excitations. The mass dependence of the phonon 
dispersion has been studied for other materials. For instance in ISj, where suppressed boron 
IE on TV of the MgB2 is related to giant anharmonicity and nonlinear electron-phonon cou- 
pling 191]. In cuprates anharmonicity could be relevant for apical oxy gen 20[ |. but for the 
in plane modes the harmonic approximation was shown to be accurate 21j and is therefore 
adopted for the present model. Moreover, as discussed in 3|, the width of the half-breathing 
phonon in cuprates is mainly due to the electron-phonon coupling. In more conventional 
materials, like a — Sn [23] or Ge 24], the isotopic frequency shifts show a dependence in- 
versely proportional to the square root of the mass and a linewidth inversely proportional 
to the mass. The error of such measurements is within a fraction of cm^^, the isotopic effect 
on the linewidth is of the order of a few cm~^ and is due to anharmonicity and isotopic 
disorder. Since in cuprates the linewidth of the anomalous phonon is an order of magnitude 
larger, one might expect that also the IE should be enhanced respectively. However, it will 
be shown below that the answer is not so straightforward. 

The quantity C,q in ([3]) contains the coupling constant, the structure factor of the BS 
phonon and the spectral weight of the response function (j2]). It can be trivially eliminated 
to obtain the ffist equation of the required set: 
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- {2^1 - 272 - ut) ■ 



(4) 
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Eq.(jl]) gives a simple approximate relation between softening and linewidth, generally a 
relatively small quantity: 7, ~ Cq {ojq — ujq) . It is useful when analyzing phonon spectra 
and making analytic estimates below. The other two equations can be obtained from ([3]) by 
calculating the isotope coefficients (IC) according to their usual definitions 

_ dlnuq _ _ dhiTq 

dlnM' dlnM' ^> 

We then get simple but somehow lengthy expressions for the two momentum dependent 
quantities (below the momentum index is omitted for brevity) that have to be obtained 
from the proposed experiments: 

a = a-0.5, h = (3-l] (6) 

V2U12 + V1U22 , V2U11 — V1U21 
a = , 



U21U12 + U21U12 + M22M11 

where 

\2 
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[g-2yf{i^^-x^ + y^f + {2xyf 

, 2y {2yf {I - gy + x^ + y^ 



{g-2y) ^g-2y)[{'^^ -x'^ + y'^f +(2yxy) 
2y 8 {yxf (1 - g^/2 + gy - x^ - y') {2yf (l - gy + + y^ 



(g - 2y) (g - 2yf ((^2 _ ^,2 + y2f + ^2yxf) (g - 2y) ((^2 _ 3,2 + y2f + ^2yxf) ' 



l-2y/{g-2y) , {2yY 2{\ ^ gy - x' - y") - g" 



(2^^) .,2 " J , :2" + 



^2 _ a;2 + y2 _ ((^2 _ 3,2 + y2f + (2|/X)') 



^2 _ ^2 + ^2 _ 2y) (^2 _ a;2 + ^2)2 + ^^yxf 



^ {2yf (l-^^ + a;2 + ^2) 3t/2 (2^)^x2 2 (1 + (71/ - ^2 - t/2) _ ^2 



- 2|/) (^2 _ 3,2 + ^2)2 ^ (2yx)' ^^2 - X2 + 2/2 _ 2y)' (^2 _ 3,2 + ^2)2 ^ 

1 2?/x2 

((7 — 2|/) ■02 — a;2 + 1/2 

All the quantities above have been scaled with fig ( Eq. (jlj) can be scaled in the same way) : 

= V<1 = TT^ ^1 = 7r^ 9g= ^- (7) 

The coefficients and 6g define the deviations from the reference values which are observed 
in many materials where such IE has been measured. It is easy to see from ([3]) that for a 
large energy separation between the phonon and electron excitations (^2^ ^ a;^) one obtains 



q;° = 0.5 and jS^ = 1 given that Uq ~ M~^/^ : these are the "standard" values [25| for the 
isotope coefficients (IC). To be mentioned that there is also a purely phononic source of the 
"normal" value (3q, since the linewidth of a phonon in isotopically pu re p erfect crystal is 



caused by anharmonicity, which scales as 7 ~ at low temperature 23|. 

When the energy of electronic excitation approaches the phonon frequency one could 
expect these deviations should grow in absolute value, and we indeed find this trend in the 
present model. However, there is a qualitative difference between the linewidth IC bg and 
the dispersion IC, a^.To clarify this point we show the explicit solution approximated by the 
first few relevant terms when Qg and Tg are the largest parameters in the problem (actual 
calculations below were carried out for the complete equations). Consequently, the linewidth 
7g and the two coefficients are small and one can obtain an analytic estimate: 

a-y i^g'^x'^ + g + '^g (y/xf -y/x^ - g^^) , 

hc^2yg{l-g-^)+2x^{l-g^l2), (8) 

where we have again dropped the g— index. We see that vanishes with phonon coupling 
constant ^g (or with the linewidth see Eq. ([3])), while the coefficient hq remains finite and 
depends on the energy separation between the two coupled excitations. Of course, when 
^g — > 0, one can not measure the mass shift of 7^, instead this relation indicates that for 
a vanishing linewidth one should expect a saturation of the respective IC at a finite value 
determined by the ratios like Xq and gq, while the dispersion IC is at its "standard" value. 
Below we will see that there is another important property of the hg coefficient, that it can 
take relatively large values even when the two interacting excitations are far from resonance. 
Here we note that the sign of the linewidth-derived IC depends on the ratio between Vtg and 

The set of three equations ((jl]) and ([5]) or the approximation (jl]) and ([8])) on the two 
parameters, Tg and is overcomplete. It therefore gives a possibility to check also the 
consistency of the model when the respective data become available. We now solve these 
equations for a, and hq with Xg and yg as independent variables to demonstrate the correlation 
between the two sets, but bearing in mind that uog and 7g are in fact intrinsically connected 
through a microscopic mechanism of electron-phonon interaction and it is the IC which are 
measured in experiments to localize the "physical" point in (xg, yg). We use the data obtained 
from inelastic neutron and X-ray scattering for some representative values of the momentum 



transfer q in Lai,s5SrQ,i^CuOi (e.g. [2, Isj) and calculate these coefficients together with 
the corresponding ratio Qq for several assumed values of the electronic frequency Vtq. At 
the F— point the frequency of the BS phonon is ujq = 85.5 meV. For momentum 1) qi = 
(0.25, 0, 0) we have Ugi = 73 meV and jgi = 13 meV; 2) for q2 = (0.22, 0, 0) we have Ug2 = 78 
meV and 7^2 = 9 meV; and 3) for ^3 = (0.2, 0, 0) respectively = 81 meV and 7^3 = 6.5 
meV. The results are shown in the table. 
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The main trends are also illustrated by the Figs. 1 and 2. As mentioned above, not all 




FIG. 1: Linewidth isotope coefficient bq calculated for = 400 meV in a broader range of x 
(= LOq/Viq) and y (= 79/^1,) . 

the values in the (x, y) plane are physically relevant and we have restricted the interval to 
the neighborhood of the data used in the table. The dispersion IC is relatively small: it 
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varies within a few percents around the "standard" value a = 0.5 even when Qq is close to 
the phonon frequency. Maximum of is reached at g = 0.25, the momentum corresponding 
to the largest softening and damping. In contrast, the IE on the phonon linewidth (bg) can 
be large even when the phonon and electron energy scales are rather distant ( 20% and 
more for the "MIR scenario"), provided the electron excitation is overdamped {g > 1). In 
this case bg is negative and ag is positive. With "standard" f3q = 1 this means a decrease 
below jSg = 0.8: as mentioned in relation to Eq.Q this coefficient tends to be negative 
when Tq > Qg and positive otherwise. Interestingly, that the sign of seems to have an 




FIG. 2: Linewidth to frequency ratio Tq/0,g of the electronic oscillator ([2]) associated to the calcu- 
lation in Fig. 1. 

opposite tendency. For a given frequency Qg the largest values of the IC bg are achieved in 
the region with maximal slope of ( 7 ) and not where softening (linewidth) has reached 
its highest value. This is because unlike a^, which mainly depends on the value of linewidth 
(proportional to softening), bq also strongly depends on the frequency ratio Uq/Qq. 

From the above analysis it can be concluded that both isotope coefficients of the half- 
breathing phonon mode contain a reach and complementary information on the electron 
excitation responsible for the phonon anomalies. This information can be used to extract 
the momentum resolved dynamical electron density response function after a detailed exper- 
imental data on the isotope coefficients become available. The response function is expected 
to be strongly inhomogeneous and anisotropic following the strong momentum dependence 
of the BS phonon softening. In our model it is fully described by the three momentum de- 
pendent parameters of the damped oscillator ([2]). If the value of bq is found to be negative (or 
/3g < 1), that would be a strong indication in favor of the "MIR scenario" and would allow to 
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study the finite momentum counterpart of the anomaly known from optical experiments by 
using the BS phonon as a probe. The accuracy of neutron and X-ray scattering experiments 
is sufficient to quantify the oxygen isotope effect, although it is "usually" weak. We have, 
however, shown that for high-Tc cupratcs the IE effect on the linewidth of the BS phonon 
mode can be large and is highly sensitive to the parameters of the electron spectrum. 

Note added by S.C. and R.C. With our deepest regret we must report that Prof. 
Maria Marinaro passed away after this work has been completed. 
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